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Abstract 

A new algorithm for extracting de Haas - van Alphen frequencies, effective masses, 
and electronic density of states contributions from calculated band energies is pre- 
sented. The algorithm creates an interpolated /c-space "super cell", which is broken 
into slices perpendicular to the desired magnetic field direction. Fermi surface orbits 
are located within each slice, and de Haas - van Alphen frequencies, effective masses 
and density of states contributions are calculated. Orbits are then matched across 
slices, and extremal orbits determined. This technique has been successful in locat- 
ing extremal orbits not previously noticed in the complicated topology of existing 
UPt3 band-structure data; these new orbits agree with experimental de Haas - van 
Alphen measurements on this material, and solidify the case for a fully-itinerant 
model of UPt3. 
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1 Introduction 



The Fermi surface of a metal is a surface of constant energy in momentum 
space (fc-space) that, at absolute zero, separates occupied and unoccupied elec- 
tron states. Determination of the Fermi surface topology of any new metallic 
compound is an important step toward understanding the physics of the mate- 
rial, since electronic properties depend sensitively on the Fermi surface shape. 
A powerful experimental Fermi surface probe is measurement of the quan- 
tum oscillatory magnetization, known as the de Haas - van Alphen (dHvA) 
effect pQ. 
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At zero temperature, in the presence of a magnetic field, the free electrons 
of a metal undergo cyclotron motion in real space. This constrains the k- 
space electron states to lie on concentric tubes, called "Landau tubes" , aligned 
along the magnetic field direction. Only the segments of tube within the Fermi 
surface contain occupied states (see Fig. [p. The radius of each tube depends 
on the strength of the magnetic field, such that as the field is increased, the 
tubes get larger and successively exit the Fermi surface. As each tube leaves the 
surface, it rapidly depopulates, causing quantum oscillations in magnetization 
and other properties. The quantum oscillations are periodic in l/B, where B 
is the magnetic field strength, with an oscillation frequency proportional to 
the extremal Fermi surface cross- sectional area perpendicular to the magnetic 
field direction pQ: 

F = (h/2ne)A (1) 



where F is the dHvA frequency, and A is the extremal area. Furthermore, as 
temperature is increased, the Fermi surface increasingly blurs, damping out 
the quantum oscillations. The damping strength depends on the effective mass 
of the electrons averaged around the extremal orbit being measured. Thus, 
dHvA measurements performed as a function of temperature and magnetic 
field direction directly reveal both the topology of the Fermi surface and the 
effective electron masses on the extremal orbits of each Fermi surface sheet pp. 




Fig. 1. Landau tubes intersecting a spherical Fermi surface. Only the sections of tube 
inside the Fermi surface have occupied states. In this case, only the orbit around 
the equator is extremal, so only one dHvA frequency would be observed. 

1.1 Basic concepts of Fermiology (may be skipped by readers familiar with 
the subject) 



For readers new to Fermi surface measurements, a few definitions are required. 
Since, at zero temperature, the Fermi surface separates occupied from unoc- 
cupied /c-space states, an orbit refers to a path along this surface that an 
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electron may trace out under the influence of a magnetic field; in our case, we 
are dealing with orbits that not only lie on the constant-energy Fermi surface, 
but on a plane perpendicular to the applied magnetic field as well. A closed 
orbit is one that forms a closed loop around some part of the Fermi surface. 
An extremal orbit is a particular closed orbit whose cross- sectional area is 
either locally maximum or locally minimum, compared to adjacent orbits on 
the same Fermi surface sheet at the same magnetic field angle. Only extremal 
orbits are detected as quantum oscillation frequencies in a dHvA experiment; 
/, II and III in Fig. [2] are examples of extremal orbits that occur at various 
angles in UPt3. 

The Brillouin zone (or equivalently, the reciprocal unit cell, depending on the 
particular volume of fc-space enclosed) is the basic unit of momentum space, 
which repeats in all directions just as the real-space unit cell of a crystal does. 
As shown in Fig. sometimes a Fermi surface links up with copies of itself in 
neighbouring Brillouin zones, forming one infinitely large, complicated shape. 
On such a surface, it is possible to have open orbits: orbits that continue 
forever in one direction, never coming back on themselves to form a closed 
loop. IV in Fig. [2] is an example of an open orbit. While open orbits cannot be 
detected in quantum oscillation measurements, their presence can be inferred 
from other experiments, such as angle-resolved magnetoresistance. Near-open 
orbits are closed, extremal orbits that would become open at a slightly different 
magnetic field angle. Ill in Fig. [2] is a near-open orbit, because a small change 
in magnetic field angle will transform it into an open orbit similar to IV. 
Finally, when determining the orbit type of a particular extremal orbit, an 
electron orbit is one that encloses occupied states, separating them from the 
unoccupied states outside the orbit; conversely, a hole orbit is one that encloses 
unoccupied states, and is associated with the motion of a "hole" through a 
Fermi sea of electrons rather than motion of a single electron itself. 



1.2 Comparison to band structure 

In order to understand the physical implications of dHvA data, the measured 
frequencies and masses are compared to those predicted by electronic structure 
calculations. However, since real compounds often possess complicated Fermi 
surfaces, including open/near-open, nested and non-central orbits, the task of 
extracting predicted dHvA orbits from calculated band energies is non-trivial. 

One previous approach involved fitting the calculated band energies to a 
Fourier series of lattice-specific star functions, which was then evaluated at 
various points in A;-space. An orbit centre and plane would be manually spec- 
ified, then the algorithm would perform a series of "stepping" and "return-to- 
surface" operations as a function of rotation angle around the orbit centre to 
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Fig. 2. The band 2 Fermi surface of UPt.3, tiled in several Brillouin zones (modified 
from [2]). I is a simple closed orbit, easily identified as extremal by visual inspection; 
II is a less-obvious extremal orbit that crosses Brillouin zone boundaries; III is a 
near-open extremal orbit; and IV is an open orbit. 

determine the cross-sectional area via Simpson's-rule integration [3]. This ap- 
proach was later expanded to included trapezoidal-rule integration for orbits 
which are multivalued at certain angles around the orbit centre [I] . While such 
an approach works well for when one knows which orbits are extremal (and 
the location of the associated orbit centres), topologically-complicated Fermi 
surfaces can have extremal orbits that are not obvious from visual examination 
of the surfaces. 



2 Numerical method 

2. 1 Overview 

Our algorithm is designed to exploit the large processing capabilities of current 
desktop computers in order to automatically extract extremal orbits, effective 
masses and density of states contributions from calculated Fermi surfaces of 
arbitrary topology, without requiring manual guidance or bias from the user. 
The code is written in the Fortran 90 language, and reads files defined in the 
XCrysDen Band-XCrysDen-Structure-File (BXSF) format. BXSF files specify 
band energies on a three-dimensional grid within a parallelepiped Reciprocal 
Unit Cell (RUC) defined by three reciprocal lattice vectors: a, b, and c [5]. 
Typical input files prepared for our algorithm contain on the order of 20 000 
/c-points, whose band energies have been calculated by an electronic structure 
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program such as WIEN2k [6J. 

Upon reading the input file, the following steps are performed: 

(1) A cubic fc-space Super Cell (SC), considerably larger than the original 
reciprocal unit cell and aligned with the desired magnetic field vector, is 
constructed (section l2.2p . A coordinate transformation maps the super- 
cell fc-point grid back into the reciprocal unit cell. Band energies at each 
of the super-cell grid points are determined from those provided in the 
reciprocal unit cell via Lagrange interpolating polynomials. The density 
of /c-points in the super-cell grid is typically much greater than that of 
the reciprocal-unit-cell grid. 

(2) The super-cell grid is divided into slices 1 fc-point thick, perpendicular to 
the magnetic field vector. On each slice, the program scans through the 
/c-points, locating as orbits the Fermi surface outlines (section I2~3"j) . The 
cross-sectional area, effective mass and orbit type (hole or electron) are 
calculated for each orbit (section 12.41) . 

(3) Orbits are matched from slice to slice, so that each orbit is associated 
with a particular Fermi surface sheet (section 12.51) . 

(4) On each Fermi surface sheet, the orbits which are extremal are singled out 
(section 12.61) . The orbit data for similar orbits found on separate sheets 
are averaged, and the results output to a file. 

(5) If automatic rotation is enabled, steps 1-4 are repeated for each new 
magnetic field vector. 

(6) The electronic density of states contribution for the band is calculated 
by comparing the reciprocal-unit-cell occupied states A;-space volume at 
energy iso-levels slightly above and below the Fermi energy to those at 
the actual Fermi energy (section 12.71) . 

2.2 k-space super cell construction 

Our algorithm operates in a large cubic fc-space Super Cell (SC), which is 
aligned with the magnetic field direction and contains many tiled copies of 
the Reciprocal Unit Cell, so the first task is to construct this cell (Fig. [3D and 
populate it with band energies. 

The reciprocal lattice vectors are defined in the BXSF file relative to fc-space 
axes xruc, Vruc, zruc, such that, for example: 

a = (a x , a y , a z ) = a x xpajc + a y VRUc + °>zZruc (2) 

According to convention, these axes are chosen so that a lies along xruc (that 
is, a y = a z = 0); for a cubic or tetragonal reciprocal unit cell, b lies along ynuc 
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Fig. 3. An example super cell, aligned with the magnetic field H, and drawn to-scale 
with the small reciprocal unit cell contained within. The super cell is cut into slices 
perpendicular to the magnetic field, which are populated with grid points. A typical 
super cell would contain 600 slices, each holding 600 x 600 grid points: far more than 
shown here. The axes for the reciprocal-unit-cell and super-cell coordinate systems 
are inset. 

and c lies along zruc- Since the super cell is aligned with the magnetic field 
(zsc\\H), it has a different coordinate system, and a coordinate transformation 
is needed to map the super-cell points back to the original reference frame: 
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(3) 



where s = sin0, t = cos0, u — 1 — cos0, v = sin 6, w = cos 6, 4> is the polar 
angle of the magnetic field measured from zruc down toward the xruc^Vruc 
plane, and 9 is the azimuthal angle of the magnetic field measured from xrjjc 
toward ymjc- Both coordinate systems share the same origin. 

No matter what the shape of the reciprocal unit cell, the super cell is always 
a cube with sides denned to be 4 x longer than the longest reciprocal lattice 
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vector. This way, for any magnetic field orientation, the super cell will contain 
enough tiled reciprocal unit cells to be able to track orbits which cross the 
zone boundaries. The super cell is situated so that 1/4 of the side length lies 
along — xsc while 3/4 of the side length lies along +xsci an d similarly for the 
other dimensions. The density of fc-points in the super-cell grid is typically 
100-200 x greater that of the reciprocal- unit-cell grid. 

The fc-points in the super cell are transformed to the reciprocal-unit-cell refer- 
ence frame via Eq. [31 and translated into the volume of the original reciprocal 
unit cell. However, since the calculated band energies are positioned uniformly 
on a reciprocal-unit-cell fc-point grid defined by the reciprocal lattice vectors, 
and these in turn are not necessarily at right angles to one another, a further 
transformation is necessary to allow the super-cell points to be compared to 
those read from the input file. A point p = p x xnuc + PyVmc + VzZruc maps 
to a point q = q a a + qbb + q c c in the input array as: 



(4) 



Once the super-cell Appoints have been appropriately mapped back to the 
space of the input array, their band energies may be derived from those orig- 
inally calculated by the electronic structure program. Since few re-mapped 
super-cell grid points will coincide exactly with reciprocal-unit-cell grid points, 
a series of third-order Lagrange interpolating polynomials [7J on a 4 x 4 x 4- 
point grid are used to determine each super-cell /c-point band energy. 

2.3 Fermi surface orbit detection 

Upon population with band energies, the super cell is cut into 1-fc-point-thick 
slices perpendicular to the magnetic field direction. In each super-cell slice, the 
program steps through the two-dimensional fc-point array, locating all Fermi 
surface orbit outlines (details of this process are shown in Fig.H]). Fermi surface 
points and associated energy slopes are determined around each orbit. Orbits 
which run into the super cell boundaries are ignored - this limits open orbits 
(the super cell is large enough that copies of non-open orbits will be found 
elsewhere). 

Note that during the "Record Fermi surface point between (x, y) and (x g , y g )" 
algorithm step (Fig. HJ), the Fermi surface point is not simply placed halfway 
between (x,y) and (x g ,y g ), but rather linearly interpolated between the two 
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, Start at bottom 
I left k-point (x,y) 
V =d4) 



Set unchecked(x,y) 
to FALSE 




Go to start of next 
row (l,y+l) 



Glance to the right 
K,y Q ) = (x+l,y) 



Set 

unchecked(x g ,y g ) to 
FALSE 




Record FS point 
between (x,y) and 
(x g ,y Q ); record 
energy slope 



Step to this glanced 
point (x,y) = (x y ) 



Glance at point 
around current point, 
clockwise from point 
we just glanced at 





Glance at point 
around current point, 
clockwise from point 
we just stepped from 



Go back to first grid 
point that was 
detected to be 
inside current orbit 



Fig. 4. Flow chart showing the Fermi surface orbit detection algorithm for finding 
Fermi surface outlines (orbits) within a given slice. Coordinates (x^y) refer to grid 
points in the two-dimensional array of Appoints on the slice. 

so that its position on the slice is given by: 



XFS 
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\E(x 



; g ,y g )-E F \(y9 y) 



(5) 



To further clarify the core Fermi surface detection algorithm presented in 
Fig. HI the action of this procedure on a small portion of an example slice is 
shown in Fig. [5j The algorithm has already been stepping around by the time 
it arrives at the depicted portion of the slice. The operations performed, in 
order, are: 

(1) Arrive at (2,3). £(2,3) is not <= E F . 
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Fig. 5. Small portion of an example slice. "+" grid points have energies greater than 
the Fermi energy; "— " grid points have energies less then the Fermi energy. Typical 
slices hold 600 x 600 grid points, so orbits usually contain many more points than 
the trivial example shown here. 



(2) Step right to the next point, (3,3). £7(2, 3) < E F . 

(3) Glance from (3,3) to the right at (4,3). -£7(4, 3) > Ep, so Fermi surface 
point a and the energy slope are recorded between (3,3) and (4,4). 

(4) Glance around (3,3) clockwise to (3,2). £7(3,2) > £7^, so Fermi surface 
point b and the energy slope are recorded between (3,3) and (3,2). 

(5) Glance around (3,3) clockwise to (2,3). £7(2,3) > £7^, so Fermi surface 
point c and the energy slope are recorded between (3,3) and (2,3). 

(6) Glance around (3,3) clockwise to (3,4). £7(3,4) is not > Ep, so step from 

(3.3) to (3,4). 

(7) Glance around (3,4) clockwise from (3,3) to (2,4). £7(2,4) > £7 F , so Fermi 
surface point d and the energy slope are recorded between (3,4) and (2,4). 

(8) Glance around (3,4) clockwise to (3,5). £7(3,5) > £7^, so Fermi surface 
point e and the energy slope are recorded between (3,4) and (3,5). 

(9) Glance around (3,4) clockwise to (4,4). £7(4,4) > £7^, so Fermi surface 
point / and the energy slope are recorded between (3,4) and (4,4). 

(10) Glance around (3,4) clockwise to (3,3). £7(3, 3) is not > £7^, so step from 

(3.4) to (3,3). 

(11) Glance around (3,3) clockwise from (3,4) to (4,3). £7(2,4) > £7^, so an 
Fermi surface point and the energy slope are recorded between (3,3) and 
(4,3). This Fermi surface point is the same as the first one found on the 
orbit (Fermi surface point a, in step 3): this orbit is done! 
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2.4 dHvA frequency, effective mass, and orbit type calculations 



For each orbit (i.e. Fermi surface outline) found in a particular slice, the dHvA 
frequency, effective mass and orbit type are calculated. Since the frequency is 
proportional to the cross-sectional Fermi surface area (Eq. [T]), it is obtained 
from the area of the polygon formed by the Fermi surface points (e.g. points 
a-f in Fig. [5]): 



h 

2vre 



N-l 



E i X FS,i UFS,i+l — XFS,i+l VFS,i) 



(6) 



where the sum is over the points on the orbit, with the last point on the or- 
bit same as the first: (xFs,N,yFS,N) = (x F s,i,yFS,i)- Note that although the 
"frequencies" are calculated for every orbit, cross-sectional areas only pro- 
duce experimentally-measurable dHvA frequencies when they are extremal; 
the extremal orbits are singled out later (section 12.61) . 



The effective mass calculation draws upon the energy slopes at the Fermi sur- 
face. As one averages around an Fermi surface outline, two different geometric 
cases are encountered. If the glance directions for point i and point i + 1 are 
parallel (e.g. points c and d in Fig. [5]), then: 
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where jfjyj is the energy slope parallel to the glance direction at point i (i.e. 

across the Fermi surface), and (j^j-) is the energy slope perpendicular to the 
glance direction at point i. If the glance directions for point % and point i + 1 
are perpendicular (e.g. points a and b in Fig. [5]), then: 
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The effective mass averaged around the orbit, given in units of the electron 
mass m e , is: 
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Note that Eq. [9] is only correct for extremal orbits, when the energy gradient 
lies within the slice. 

Since the orbit-finding algorithm (section 12.31) always steps around the inside 
of electron orbits (e.g. Fig. E]) and the outside of hole orbits, the orbit type is 
determined by comparing the orbit area to the area bounded by the polygon 
formed by the "stepped-on" grid points. If the orbit area is larger than the 
stepped-on area, it is flagged as an electron orbit; if the reverse is true, it is 
flagged as a hole orbit. At this point, average coordinates, coordinate standard 
deviations, maximum coordinates, and minimum coordinates for the orbit 
outline are also calculated. 



2.5 Slice-to- slice orbit matching 



Once the orbits outlines have been located on all slices in the super cell, Fermi 
surface sheets must be reconstructed by matching orbits on adjacent slices. In 
order for an orbit i on one slice to be matched with an orbit j on the preceding 
slice, and therefore added to orbit j's sheet, all of the following conditions must 
be met: 

• The average xfs and y F s coordinates of orbit i are both within one standard 
deviation of the average xfs and yps coordinates of orbit j. 

• The maximum xps an d Vfs coordinates of orbit i are both within two 
standard deviations of the maximum xps an d Vfs coordinates of orbit j. 

• The minimum xps an d Vfs coordinates of orbit % are both within two stan- 
dard deviations of the minimum xfs and Vfs coordinates of orbit j. 

If multiple orbits on one slice satisfy the conditions for matching with an orbit 
on the preceding slice, the parameter Bi is calculated for each candidate orbit, 
and the orbit with the lowest B~ value is chosen for the match: 



B: 



+ 
+ 



1 z r 

avg [XFsJj ~ avg (xFs)i\ + [avg (y F s)j ~ avg (y F s. 

i 2 r 

max (x F s)j ~ max (x F5 )J + [max (y FS ) j - max (y FS 



(10) 



min {xps)- ~ min (xps), + niin (y F s) j ~ min (y F s)i 



2.6 Extremum determination 



If an orbit has a frequency (cross-sectional area) which is greater than both of 
the adjacent orbits on the same sheet, or less than both of the adjacent orbits, 
it is taken to be extremal. Once all extremal orbits in the super cell have 
been selected, they are sorted by frequency, and those which have frequencies 



11 



within 1% of the next smallest orbit are taken to be multiple copies of the 
same orbit. Since the super cell contains more than one reciprocal unit cell, 
it is expected that multiple copies of a given extremal orbit will be found. 
The copies have their frequencies, effective masses and orbit types (electron 
= 1, hole = -1; useful for confirming that orbits of different type haven't been 
put on the same sheet) averaged, with standard deviations used to provide an 
indication of error bounds on these quantities. 



2. 7 Density of states calculation 



Calculation of the band contribution to the electronic density of states pro- 
ceeds separately from the super cell calculations, and is only performed once 
for a given input file (compared to the many repeated super-cell iterations for 
different magnetic field directions). To determine the density of states contri- 
bution, the following steps are executed: 

(1) The band energies within the reciprocal unit cell are interpolated very 
finely. No super cell is used. 

(2) The /c-space volume of the reciprocal unit cell and volume per interpolated 
/c-point are calculated. 

(3) The number of occupied points (those with < Ep) are counted, and 
the occupied volume Vp calculated. 

(4) The number of points within an iso-surface at E + = E F + dE, with 
typical dE = 0.00001 E F , are counted, and the volume of such points V + 
calculated. 

(5) The number of points within an iso-surface at EL = Ep — dE are counted, 
and the volume of such points V- calculated. 

The electronic density of states (DOS) contributions are then: 



\V F - V+\ 



with an indication of error bounds given by the difference between + and 
values. 
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3 Results 



3.1 Test Fermi surfaces: sphere and cylinder 



In order to validate our algorithm, two artificial cases with exactly known 
results were constructed and used as program inputs. First, a spherical Fermi 
surface was generated, with band energies that depend on the square of the k- 
space distance from the Fermi surface centre. The input file contained band en- 
ergies specified on a 28 x 28 x 28 grid in a cubic reciprocal unit cell, picked such 
that the effective mass m* = 10.543 m e and dHvA frequency F = 1.568 kT, 
constant for all magnetic field angles. Second, a cylindrical Fermi surface was 
generated, with band energies that depend on the square of the fc-space dis- 
tance from the zruc axis. The input file energies were specified on the same 
grid and reciprocal unit cell as the spherical case, and picked such that the 
effective mass m* = 3.457 m e and dHvA frequency F = 2.346 kT when the 
magnetic field is aligned with the Zruc axis. As the field is rotated away 
from the i^uc axis by an angle (ft, the frequency should increase as 1/ cos</>, 
asymptoting to an open orbit as (f> — > 90 degrees. 




Fig. 6. de Haas - van Alphen frequency plotted as a function of magnetic field angle 
4> for spherical (bottom) and cylindrical (top) test Fermi surfaces. The cylinder axis 
lies along the (ft = direction. Data extracted by our algorithm are shown as points; 
exact results are shown as solid lines. The inset shows a close-up of the low- angle 
region, as well as the exact input effective masses and corresponding band masses 
extracted by our algorithm for both cases. 
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Fig. [6] shows the results our program obtained from these two test cases. In 
both greement with the input values is excellent, particularly in light 

of the coarseness of the starting band energy grid. For the spherical case, our 
algorithm extracted a frequency of F — 1.566 kT and mass of m* = 10.539 m e , 
varying by less than 0.002 kT and 0.002 m e respectively over the angular range 
tested. For the cylindrical case, with field directed along the zrjjc axis, our 
algorithm extracted a frequency of F = 2.345 kT and mass of m* = 3.460 m e . 
As the field is rotated away from the z^juc axis, the extracted frequency follows 
the predicted 1 / cos 4> curve very closely. 



3.2 UPt 3 



Having validated our computational approach, we applied it to a material of 
current interest: the heavy fermion superconductor UPt 3 , which has a hexago- 
nal crystal structure. The term "heavy fermion" is used to describe this mate- 
rial, because strong electron-electron correlations lead to a very large effective 
mass - the conductance electrons act very "heavy". Traditional UPt3 band 
structure calculations, in which all uranium 5 /-electrons are itinerant rather 
than localized, find 5 energy bands crossing the Fermi energy [8], including 
the complicated Fermi surface sheet shown in Fig. [2j Comprehensive quantum 
oscillation experiments have been performed on this compound [2], and while 
the orbit shapes generally agree well with band structure, some frequencies 
seen experimentally had no counterparts from the calculation. 

By using the previously-calculated UPt 3 band structure information [5J as in- 
put to our program, we were able to rigorously determine the frequencies that 
should be seen by dHvA. Through this process, in addition to confirming the 
previously-known orbits, we found new extremal orbits in the old band struc- 
ture calculation that correspond to the experimentally-measured frequencies 
which had until now been orphaned [2]. One such frequency, labeled rj, corre- 
sponds to orbit II in Fig. [2]- since it is difficult to see by eye that this orbit 
is extremal, it had not been noticed prior to our work. 

Full frequency versus magnetic- field- angle results for all 5 -Ep-crossing bands 
are shown in Fig. [7J Good agreement between measured and calculated Fermi- 
surface shape (manifested in the frequency variation with angle), but not size 
(manifested in the overall magnitude of the frequencies) is usual for heavy 
fermion materials, because their unusually flat bands mean that small shifts 
of the Fermi energy can produce large changes in the size of a Fermi sur- 
face sheet. Since our calculation has located previously- missing orbits in the 
traditional band structure data, thus bringing it into closer agreement with 
the experimental results, it has helped resolve a recent controversy over the 
fundamental itinerant nature of the UPt3 5/ conduction electrons [2]. 
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Fig. 7. Measured and calculated UPt3 de Haas - van Alphen frequencies as a function 
of magnetic field angle (from |2j). Solid and dotted lines are frequencies extracted 
from a fully-itinerant band structure calculation by our algorithm. Points are ex- 
perimental measurements, with Greek letters indicating orbit labels. The r] orbit 
(orbit in Fig. [2j) is one of the new extremal orbits found in the band structure 
data by our algorithm. 



4 Conclusion 



We have developed a new approach for extracting quantum oscillation fre- 
quency, effective mass, orbit type and electronic density of states contribution 
from calculated band energies. By employing a large, heavily-interpolated k- 
space super-cell and exploiting recent advances in desktop computing power, 
our program can robustly characterize complicated Fermi surfaces to locate 
all extremal orbits. 

In addition to excellent performance on simple test Fermi surfaces, when ap- 
plied to the complicated energy bands of UPt 3 our algorithm located new 
frequencies which had not been previously found on the calculated Fermi sur- 
face. These new frequencies increase the agreement between model and exper- 
iment, providing new evidence for the itinerant nature of the 5/ electrons in 
this material. 
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